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Abstract. We consider the gravitational potential and the gravitational rotation field generated by an spherical mass 
distribution with exponential density, when the force between any two mass elements is not the usual Newtonian one, 
but some general central force. We invert the usual integral relations and obtain the elemental interaction (between 
I two point-like masses) as a function of the macroscopic gravitational field (the one generated by the distribution). 
■ Thus, we have a direct way for testing the possibility of finding a correction to the Newtonian law of gravity that can 
^j" explain the observed dynamics at large scales without the need of dark matter. We show that this method can be used 
even in the case of spiral galaxies with a good level of confidence. 
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Q 1. Introduction 

h ; 

The dynamical analysis of rotation curves of galaxies, binary galaxies, clusters of galaxies and the structures known 
C$ at large scales show large discrepancies between the observed behaviour and the one expected from the application 

of General Relativity and its Newtonian limit to the visible mass. This disagreement has led many astrophysicists to 
• »"H believe in the existence of a large amount of non-visible matter and is, thus, commonly known as the Dark matter 
?\ problem. 

In spite of this, there is no direct evidence for the validity of either Einstein's General Relativity or Newtonian 



gravity at scales much larger than those of the Solar system (See, e.g., Will 1993). There is therefore no experimental 
or observational reason to ascertain that unmodified General Relativity holds at larger distances. This leads us to 
think that we should be open to the possibility that it had to be revised (perhaps in the same spirit as Newton's law 
had to be modified for strong fields and large velocities). 

In this paper we consider the possibility that Newton's law of gravity is just a good approximation at short distances 
of a more general expression for the force. It is interesting to identify which, if any, extensions of the usual inverse 
square law are compatible with the dynamics observed at large scales. 



Work along these lines has already been done (Tolhine 1983 , Kuhn & Kruglyak 1987 , Mannheim & Kazanas |l989 ) 



assuming a specific functional form for the force, and then evaluating the field generated by a mass distribution (for 
instance, a galaxy) by performing the corresponding three-dimensional integrals. We present and work out a method 
that allows us to follow the inverse methodology, that is, to infer, directly from observations, the phenomenological 
law of gravity that is able to generate a given macroscopic gravitational field. We do it for the case of a spherical mass 
distribution with exponential density. 

In Section || we give the general definitions that will be used later in Section || for the case of spherical symmetry with 
exponential density. In Section [| we show the mathematical basis underlying the results presented in Section ||. Finally 
we offer some conclusions. In Appendix ^ we consider the possibility of applying these results as an approximation 
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to the study of spiral galaxies overcoming the fact that they do not show spherical symmetry. In Appendix |b] we list 
some of the mathematical identities used in Section ^|. 

2. General definitions 

Let us assume that the gravitational potential generated by a point-like mass does not correspond to the usual 
Newtonian form but can be written in terms of a function g(r) that describes the deviation from the Newtonian law, 
that is, 

_ G mim 2 

<j>(r) = g{r). (1) 

r 

where <p(r) is the gravitational potential experienced by two point-like particles of masses mi and mi separated by a 
distance r and Go is the Newton's constant. 

This modification could, for example, be due to the many body nature of the mass distribution making up the 
galaxy, quantum corrections, a relativistic theory different from General Relativity... This is irrelevant in what follows. 

The force per unit mass is, by definition, the gradient of the potential, 

G 771lTO 2 r 

F(r) = - g eS (r)-, (2) 

where we have introduced 

9cs(r) = g(r) - rg'(r). (3) 

In this way, to find the total potential or the total force generated by a mass distribution f2 with density p(r), one 
must integrate over the volume spanned by Q to get: 

for the potential experienced by a point mass at a distance R from the centre of f2, and 
for the force. 

In the case that the gravitational potential is only a function of the distance to the centre of the distribution, it is 
convenient to introduce two new functions ^(R) and 'I'eff(-R) such that: 

m - - *(*), w 



F ( R ) = *off(i?) ^, (7) 

and the rotation velocity of a test particle in a circular orbit bound to the distribution will be: 

O*) - ^f* *e ff (i?) (8) 

where the auxiliary functions ^off (R) and ^(R) satisfy the following functional relationship: 

* eS (R) = ¥(R) - RV(R), (9) 

Our goal is to design a procedure where, assuming that F(R) is known (say from observation of the rotation 
velocity) for all values of i?, we obtain a g e s(f) that generates the given rotation velocity. Or, in other words, given 
the potential as inferred from observations we want to find which g(r) could have generated it. Actually, what we will 
find are g(r) and g c s(r) as functions of ^ c s(R) and *&(R) respectively. 
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3. Spherical mass distribution with exponential density 

In view of the applications that we have in mind, we study a spherically symmetric distribution with an exponentially 
decaying density: 

p(r) ee Po e- ar . (10) 

Our ultimate goal is to find a method to study the discrepancies between the observed rotation curves of spiral 
galaxies and the curves predicted by using Newton's law of gravity. The luminosity profile of many spiral galaxies can 
be well fitted ass uming that the density of luminous matter decreases exponentially with distance from the centre of 
the galaxy (Kent |l987|) . This is the reason why we are interested in studing such a density function, even though spiral 
galaxies are not spherical. 

Using Eqs. (|), (|) and ©, in Eqs. (@) and |), and considering spherical symmetry, the two problems sketched 
in Section ^| can be conveniently recast in the form of two integral equations: 

(i) Given ^(R), find a function g(r) that satisfies the equation: 

r* r*t r^ RA+r ''- 2Rrc ^ si nee-- = %m (11 ) 

Jo Jo Jo VR 2 +r 2 - 2Rrcos8 a 3 R 

and 

(ii) Given ^ Q g(R), find a function g e s(r) such that: 

dr r + T% - 2B ^1(R - rcoseysinOe- = (12 ) 

o Jo Jo (R 2 + r 2 - 2Rrcos8)i a R 

In the next section the solution to these integral equations will be described in detail. The results can be summarized 



(i) Potential problem (viz. Eqs. (|6J) and (|l l[) ) 

In this case, the exact solution to the problem is 
g(x) ^"(x) + (x) (13) 

where the function \1/ has the following behaviour at the origin: 

*(0) = *"(0) = 0. (14) 

(ii) Force and velocity problem (viz. Eqs. (0), (||) and fll^)). 

Here, the exact solution is given by the following expression: 

9cS ( x ) = 0, eS ( x ) - ±*^(x) + ±9%\x) + - -±- A [2x% s (x) - 2x 2 ^ s (x) + x 3 ^(x)\ (15) 

or a 4 orx a 4 a; 4 

The behaviour of \& at the origin is as follows: 

* cff (0) = *^(0) = < ff (0) = 0. (16) 

The behaviours at the origin just tell us that ^ c g(R) oc R 3 for R ~ 0, and thus, V ro t(R) oc R for R ~ 0, which is 
in fact in good agreement with the observations (as the observed rotation curves are usually well fitted in the inner 
regions by a straight line) and, for a non-Newtonian gravity point of view, it is also in agreement with the fundamental 
experimental constrain that for short distances the gravitational interaction must be well described by a Newtonian 
limit. 

4. Mathematical development 



First, we will show the solution to what we call the potential problem, that is, how to go from Eq. (|ll|) to Eqs. (13) 



and (|14|). Later we will use these results to solve the force problem, i.e, how to go from Eq. ( |12|) to Eqs. ( pL 5| ) and ( pL 6[ ) 
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4-1- The potential problem. 

In order to go from equation Eq. ( |Tl"| ) to Eqs. ( |l3| ) and (14), we will need to use several mathematical identities. For 
convenience, these are listed in Appendix M. 



Inserting Eqs. (Bl), (E33h, (B5), (B6) and (B7) into Eq. (O) we can write as: 



dp- 



M 
p 



dre 



■ u a3 d 
r smpr smpR = — 

7r da 



, 9 sip) . D 
dp^ o sin pit 

TVS I /-..i 



or 



(17) 



Then, we can apply Eq. (B2) and invert the Fourier transform to obtain, after some straightforward calculations, 
the following more useful form: 



*(i?) = -f {(1 - aR)e aR dx g{x) e~ ax - / °° dx g[x) e" tta 



(l + aR)e- aR J Q dx g(x) e ax - dx g{x) e~ ax 



+ 



(18) 



a aR 



dx g(x) x e 



J °° dx g(x) x e 



ae 



-aR 



dx g(x) x e ax + J °° dx g{x) x e ax \ 



In order to simplify these expressions, we introduce an auxiliary function ip(x) that makes the integrals exact: 



g(x) = ip(x) - \ ip"ix) + \ ip m ix) 



a* 



a* 



We can insert Eq (|n]) into Eq. (Q) and, upon integration by parts, we get: 



*(iZ) = - (1 + ^-)e~ aR m + A 



^"(0) 



(19) 



(20) 



where ^ is a solution to the ordinary differential Eq. (|l^) that satisfies the conditions of being an analytic function at 
x = 0, and 

lim V (fe) (z) a exp(-aa;) =0 ; fc = 0,l,2,3 (21) 

x — >oo 

where ip( k \x) stands for ip{x) and its three first derivatives. 

These conditions are easily fulfilled in all the cases of interest. Actually, the analiticity is satisfied in the Newtonian 
limit, that is the behaviour that we expect to recover at R ~ 0. Although it is possible to artificially build a pathological 
ipix) such that it can represent a physical system without satisfying Eq. (pl|), it can be seen that almost every tp function 
that does not satisfy it corresponds to a rotation velocity that grows almost exponentially with the distance, which 
clearly seems to contradict the observations. 

Actually, it is straightforward to see that, provided ip{R) is a solution to Eq. ( ]19|) , then ^(i?) is also a solution 
to the same equation. Moreover, the terms proportional to V'(O) and ^"(0) in Eq. ( P0[ ) assure that ^> and its second 
derivative are both zero at the origin. Taking all of this into consideration, we finally obtain that: 



g(x) = *(x)- ^U"(z) + ^ iv Hx), 



*(0) = W"(0) 



(22) 
(23) 



4-2. The force and the rotation velocity. 

Once ^(-R) is known we can calculate ^> e giR) using Eq. @. Equivalently, once g(r) is known, g e sir) can be obtained 
through Eq. Using these two equations together with Eq. (^|]), and after some straightforward calculations, we 
can find a direct relation between g e g and ^cff^ 



g e six) = ^ e s(x) 



-[2x% H (x) - 2x 2 ^ s (x) + x^ix)) 



(24) 



■ , II >•• / j -' rl 2 ^esi x ) 

a a ax u-ju- 

From the behaviour of \& at the origin (Eq. (p3|), and Eq. (||)), it is easy to see that, at the origin, will satisfy: 

* cff (0) = <ff(0) = < ff (0) = 0. (25) 
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5. Discusion and conclusions 

We have found the exact solution to the problem of inverting the integral relation between the elemental law of gravity 
and the gravitational field generated by a spherical mass distribution. We have arrived at a direct and simple expression 
which can be useful in order to infer the law of gravity that could explain the large scale gravitational behaviour if 
we had good data on rotation curves of spherical galaxies with an exponential density profile. However, most of the 
observational data available on rotation curves are for spiral galaxies. Despite the fact, shown in Appendix |X[ that this 
formalism is a good approximation to study this class of galaxies if g c s is a growing function of r, we believe that it is 
better to study the case of a thin disk distribution that is a much better approximation to the real morphology of spiral 
galaxies. In two forthcoming publications we will give a similar solution for a thin-disk distribution (Rodrigo-Blanco 



1996) and we will apply it to real galaxies (Rodrigo-Blanco & Perez- Mercadcr 1996). 



A. Appendix: Gravitational difference between a disk and a sphere 



|i = 0.0 (Newt) \L = 0.5 n = 1.0 




0.00 1 — ^ — — ^ — 1 0.0 1 ■ 1 ■ 1 — ^ — 1 ■ 1 0.0 1 ■ 1 ■ 1 — ^ — 1 ■ 1 

02468 10 2 4 6 8 10 02468 10 



H = 1.5 n = 2.0 n = 3.0 




02468 10 02468 10 02468 10 

cxR ciR aR 

Fig. 1. Rotation velocities for a sphere (solid line) and for a disk (dots) when g c s(r) = (^) M for some values of ji. In every 
case, for the sake of clarity, the velocities are normalized by dividing by the appropiate constant G °,^f^ a ■ 



At first sight, we would think that we cannot apply the results of the previous sections to spiral galaxies, because they 
have been obtained assuming spherical symmetry, and spiral galaxies are disk-shaped, not spherical. In spite of this, 
it seems logical that the difference between the gravitational field generated by a sphere and the one generated by a 
disk becomes negligible at very large distances (here "large distances" means large compared with some typical length 
scale for the distribution; it could be, for instance, a -1 , i.e, the exponential length scale of the mass distribution). 
Furthermore, when we consider the effect of an increasing g c s{ r )y it is evident that the meaning of what is a "large 
distance" changes as we change the functional form of g e s(r). That is, the faster g e s{r) grows as a function of r, the 
smaller the distance at which a sphere is indistinguishable from a disk becomes, when considered from the gravitational 
point of view. 
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In order to obtain a more quantitative idea of the gravitational difference between a sphere and a disk, when g c s 
is a growing function of its argument, we have chosen a parametric family of g c g 's given by: 

g cS Jr)^ (£f (Al) 

where p parameterizes how fast c/eff grows. 

Then we have calculated the rotation velocity performing the three-dimensional integral numerically for several 
values of p. We have done it both for the case of spherical symmetry, and for the case of cylindrical symmetry. We 
have assumed for the disk a small thickness of ft, = aT 1 /6, and the rotation velocity is calculated in the plane of the 
disk. 

In figure ([!]) , we have plotted both velocities for a disk as well as for a sphere for some values of p. In each case we 
have normalized the solutions dividing by a convenient constant Vq defined as: 

It can be seen that both curves tend to merge as p increases. 

However, we would like to have a more quantitative way of describing the difference between both curves as a 
function of p. In order to do this, we define a quantity a 2 -, (see below for its meaning) as follows: 

^) = lv^ ^ (A3) 

Here the subscripts D and S stand for 'disk' and 'sphere', respectively. We sum over fj, which are the points where 
the integrals are calculated. The total number of points for each value of p is N = 100. 

Because of the way it is defined, cr|> (p) is a measure of the mean square error that we make in the rotation velocity 
if we consider a sphere instead of the real disk, for each value of p. We thus aim for a value as small as posiblc for 
crfj. In Figure (||), we plot the value of a 2 D versus p. Once again, it can be seen that the larger the value of p is, the 
smaller the difference. 

Considering what these plots mean, we see that our mathematical results can be used in the case of spiral galaxies, 
with a good level of confidence, provided g e s{f) grows fast enough with r. Actually, the most popular corrections to 



<7 e ff = 1 that can be found in the literature are g c s oc r, (see Tohline 1983 and Khun & Kruglyak 1987), and g c ^ 



oc r 



2 



(see Mannheim & Kazanas 1989). It can be seen in the figures that the approximation is quite good in both cases. 
Although the results obtained in this way are not exact, they give a qualitatively correct picture. However, it will be 
better to use the solution to the problem in the case of disk symmetry itself, as we will do in a forthcoming publication 



(Rodrigo-Blanco 1996) 



B. Appendix: Mathematical identities 

1. Fourier sine transform 



2 r°° 

g{\fr 2 + R 2 - 2rRcos9) = - / g s (p) sm{p^/r 2 + R 2 - 2rR cos 9) dp (Bl) 

f Jo 

/•OO 

9s(p) = / g(t)smptdt (B2) 



2. Addition theorem for Bessel functions (see Gradshteyn |1980[ ). 

M^) = ilrfr,) f> + k ) J ^ {mp) J - +k(mr) q;(cos e) (bs) 

k=0 1 

where: 
lo = \J r 2 + p 2 — 2rp cos 9 

p<r (B4) 
C£ = Gegenbauer Polynomials. 



Using Eq. ([f33|), together with: 
sin(mz) = J^^-Ji/ 2 {mz) 



(B5) 
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1.0 



Fig. 2. Mean square error (o^) in the velocity as a function of p, when a sphere is considered instead of a disk, and using 
g e ff(r) = The exact meaning of crfj is explained in the text. 



we get: 
sin 



pVr 2 + R 2 - 2rR cos ( 



Vr 2 + R 2 -2rRcos6 
3. Orthogonality of Legcndrc Polynomials. 

d0P k (cos9) sin (9 = 2<5 k ,o 



-j= J2 (2k + 1 ) J k+ 1 (pr ) J k+ 1 (pi?) P fc (cos 0) 
v Hr k=0 



/ 

Jo 



(B6) 



(B7) 
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